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3-D mapping of diffuse scattering in Pb(Zn!/ 3 Nb2/3)03-xPbTi03 
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High energy (67 keV) x-ray diffuse scattering measurements were performed on single crystals of 
Pb(Zn 1/ / 3 Nb 2 /3)03-a;PbTi03 (PZN-xPT). A novel technique was developed to map out the diffuse 
scattering distribution in all three dimensions around a large number of Bragg peaks simultaneously, 
taking advantage of the almost flat Ewald sphere of the high energy x-ray beam. For x = 0, 4.5%, 
and 8%, the results are very similar, indicating same type of correlations of polarizations in these 
compounds. Our results show that the diffuse scattering intensity consists of six (110) rod-type 
intensities around reciprocal lattice points. A simple model is suggested where {110} type planar 
correlations of the in-plane (110) type polarizations in real space contribute to the (110) rod-type 
diffuse intensities in the reciprocal space. The planar correlations of polarizations are likely a result 
from condensations of soft (110) polarized optic phonon modes. 



PACS numbers: 77.80.-c, 77.84.Dy, 61.10.Nz 
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I. INTRODUCTION 

Relaxors are a special class of ferroelectrics whose di- 
electric properties show diffusive and strongly frequency 
dependent phase transitions. Pb(Zn 1 / 3 Nb 2 /3)03 (PZN) 
is one of the prototype relaxors, which has attracted 
much attention due to its high piezoelectric properties 
when doped with PbTiOa (PT)iSii. One unique prop- 
erty about relaxors is the appearing of local polarized 
nano-sized regions at the Burns temperature Td, which 
is a few hundred degrees above the Curie temperature 
Tc- This was first suggested by Burns and Dacol^ in 
interpreting their measurements on the optical index of 
several relaxor systems, including PZN and one of its 
close analog, Pb(Mg 1/3 Nb 2/3 )03 (PMN). 

The PNR was then extensively studied by diffuse scat- 
tering measurements. It was found that in PMN, the 
diffuse scattering starts to appear at around Td « 600 K, 
and increases monotonically with cooling^Si^. Similar 
results were also observed in PZNS (see Fig. [TJ, where 
Td is much higher than that of PMN. In order to under- 
stand the nature of polarizations and correlations/shapes 
of the PNR in relaxors, one of the most direct methods 
is to study the spatial distribution of the diffuse scat- 
tering intensity. Recently, several neutron2*ifi*ii*iSiii and 
x-rajii^i diffuse scattering measurements have been car- 
ried out on various relaxor compounds to investigate the 
shapes of diffuse scattering intensity distributions in the 
reciprocal space around different Bragg peaks. 

However, the majority of these measurements were 
performed in a 2-D scattering plane, i.e., probing the two- 
dimensional diffuse scattering distribution around Bragg 
peaks. For example, measurements for PMN around 
(100) and (110) Bragg peaks in the (HK0) scattering 
plane (see the schematic representation in Fig. ^) sug- 
gest that the diffuse scattering intensity extend along 
the (110) directions. But there is still ambiguity about 
whether or not these are projections of diffuse scattering 
intensities along other out-of-plane directions (e.g., pro- 
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FIG. 1: Neutron diffuse scattering intensity measured around 
(100) and (300) peaks for PMN and PZN, by Xu et al& and 
Stock et aim, respectively. The small anomaly in the PZN 
diffuse data around T — Tc is likely due to critical scattering 
at the phase transition. The inset shows schematics of diffuse 
scattering intensity distributions in the (HK0) plane, around 
the (100) and (110) Bragg peaks. 



jection of [111] direction) onto the scattering plane. You 
et alm^ were the first to measure the out-of-plane x-ray 
diffuse scattering intensities of PMN around the (300) 
Bragg peak. Yet no systematic 3-D mapping of diffuse 
intensity distributions around different Bragg peaks in 
these relaxor systems has been reported up to date. The 
lack of understanding of the diffuse scattering distribu- 
tion in three-dimensions made it very hard to establish 
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an effective and convincing model that can explain and 
predict accurately the diffuse pattern. 

In this paper, we present high energy x-ray diffraction 
results, on the diffuse intensity distribution in three di- 
mensions measured simultaneously around many Bragg 
peaks of PZN-xPT single crystals. Our results show that 
the diffuse scattering intensities extend along all six (110) 
directions around most Bragg points. However, for cer- 
tain Bragg peaks, one or more of the six (110) type dif- 
fuse intensities are absent, because of the associated po- 
larization (atomic shift) being perpendicular to Q. We 
suggest a simple phenomenological model, where (110) 
polarizations are planar correlated in the {110} planes. 
For example, [110] polarization correlated in the (110) 
plane. The diffuse scattering patterns calculated based 
on this simple model are in good agreement with most 
previous and current results of measurements on diffuse 
scatterings from PNR in relaxors. Details of this novel 
technique, our results, and the model calculations, are 
discussed in the next three sections. 
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II. EXPERIMENTAL TECHNIQUES 

Single crystals of PZN, PZN-4.5PT (4.5PT) and PZN- 
8PT (8PT) have been studied. The PZN single crystal 
is 3 x 3 x 1 mm 3 is size, and was grown at the Simon 
Fraser University in Canada (the same crystal used in 
Ref[H). The 4.5PT single crystal is 5 x 5 x 3 mm 3 in 
size and provided by TRS ceramics. The 8PT crystal is a 
part (~ 4 x 4 x 1 mm 3 ) of the original crystal previously 
studied by Ohwada et alH, grown at the Pennsylvania 
State University. 

The x-ray diffraction measurements were performed 
at X17B1 bcamlinc of the National Synchrotron Light 
Source (NSLS). A monochromatic x-ray beam of 67 keV, 
with an energy-resolution of 10~ 4 (AE/E), was produced 
by a sagittal-focusing double-crystal monochromator us- 
ing silicon [311] reflection with both crystals in asymmet- 
ric Laue mod©i&. In principle, with a four-circle x-ray 
diffractometer, one can reach any points in the recipro- 
cal space, if allowed by geometry. However, in reality, 
most measurements on diffuse scatterings are performed 
in the "zone" . In other words, the measurements are per- 
formed mostly in a plane - the diffraction plane, around 
certain Bragg peaks. This is also the case for neutron 
diffuse scattering measurements. For example, measure- 
ments in the (HK0) zone indicate that the measurements 
are performed in the plane perpendicular to the c-axis (of 
the pseudocubic system), so that L is fixed (to 0). This 
only probes two of the three dimensions of the diffuse 
scattering distribution around Bragg peaks. Certainly 
one can measure diffuse scatterings in other planes and 
essentially reconstruct a 3-D picture. But the beam time 
and efforts spent would usually be enormous. 

Position sensitive detectors such as CCD detectors or 
Image-plates have been used to perform monochromatic 
Laue style measurements. In this type of measurements, 



FIG. 2: Schematics of the high energy x-ray diffuse measure- 
ment set-up. The dots denote the reciprocal lattice points, 
and the arc describes where the Ewald sphere cuts in recip- 
rocal space. The diffracted x-ray beam will then go into the 
CCD camera as shown after the sample. 



scattering intensities in a plane, i.e., the Ewald sphere, 
can be taken simultaneously, very similar to the case of 
electron diffraction measurements in transmission elec- 
tron microscopies. This usually requires Bragg points to 
be close to the Ewald sphere. In some cases, single crystal 
samples were oscillated while the image was taken. This 
is a very efficient way to obtain 2-D diffraction intensity 
distributions, but without any information, or with very 
coarse resolution, along the incident beam direction per- 
pendicular to the image plane. One example of this t ype 
of measurements on relaxor systems is shown in Ref. Il9t 
In our measurements, we were able to take advantage 
of the high energy (67 keV) of the x-ray beam and a 
position sensitive CCD detector to perform 3-D x-ray 
diffuse scattering measurements in a very efficient man- 
ner. The focusing double-crystal monochromator set-up 
is also essential since it produces more than two orders of 
magnitude more photon flux than that of a conventional 
set-up. A schematic of the experiment set-up is shown 
in Fig. |2 The attenuation length of 67 keV x-rays in 
our samples are a few hundred microns, and the mea- 
surements are performed in a transmission mode. The 
reciprocal lattice and the radius of the Ewald sphere (the 
curve in Fig. [21 are scaled from the real values in our mea- 
surements (a* = 27r/a = 2?r/4.06 A = 1.5476 A" 1 and 
h = 33.87 A" 1 ). The incident x-ray beam was assumed 
to be along the reciprocal lattice [001] direction, and the 
[100] direction is pointing up. Because of the high x-ray 
energy, and thus a large ki, the Ewald sphere is almost 
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flat, parallel to the (001) plane, perpendicular to c-axis at 
small Q. This makes it possible to reach close to many re- 
ciprocal lattice points in the (HKO) zone simultaneously 
with one sample geometry. In addition, it is important 
to note that at larger Q values, the Ewald sphere is not 
exactly cutting through the reciprocal lattice points. In- 
stead, the sphere cuts through the reciprocal lattice at 
none-zero L values, as well as a small tilt angle to the 
(001) plane. For example, around the (200) reciprocal 
point, the Ewald sphere is about SL — 0.09 r.l.u. away 
from the exact (200) position. So we are in fact measur- 
ing diffraction intensities at (H,K,-0.09) around the (200) 
position on the CCD. By tilting the sample, therefore tilt- 
ing the reciprocal lattice, the Ewald sphere can cut the 
reciprocal lattice at different SL values as we need. With 
a few measurements taken at different sample tilts, the 
geometry of the diffuse scattering intensity distribution 
in three dimensions Idiff(H,K,L) can be mapped out 
easily around selected Bragg peaks. Of course, in our 
measurements, the Ewald sphere is curved not only ver- 
tically, but also horizontally so that the situation would 
be very similar around the (020) peak as the (200) peak. 

Our measurements show that results from PZN, 4.5PT 
and 8PT crystals are qualitatively the same. The x-ray 
beam size used in the PZN and 4.5PT measurements were 
0.5 x 0.5 mm 2 , and 0.2 x 0.2 mm 2 for the 8PT measure- 
ments. The g-resolution of our measurements is mainly 
limited by the x-ray beam size. The 8PT measurements 
with finer resolution show the features in a better de- 
tail, so in the next section we will focus on the results 
obtained using the 8PT sample. 



III. 3-D X-RAY DIFFUSE SCATTERING 
MEASUREMENTS 

A. (HKO) zone 

Fig. |3 shows a CCD image taken at T = 200 K, 
which is well below the ferroelectric transition temper- 
ature T c ~ 450 K for PZN-8PT2Bi2i. This was measured 
in the (HKO) zone, as shown by the schematics in Fig. El 
The x-ray beam is incident in the [001] direction, perpen- 
dicular to the CCD plane. [100] and [010] directions are 
marked in the figure. The center of the image, which is 
the origin of the reciprocal lattice, is black because a cen- 
tral beam stop was used to block the direct beam. The 
tilt (vertical) and rotation (horizontal) of the sample has 
been aligned carefully so that the intensity around (100), 
(100) are the same, as well as (010) and (010), to ensure 
that the CCD plane is as parallel to the (001) plane as 
possible. The counting time was 200 seconds, and we 
were able to measure very fine features simultaneously 
around many Bragg peaks for Q as large as (2,2,0). The 
rings in the image are powder rings from the Be win- 
dow in the high temperature displex, and the rest of the 
background is mainly due to air scattering. 

At the four {100} positions ((100), (010), (100), and 



(010)), the Ewald sphere is still close to the Bragg peak 
(SL ~ 0.02 r.l.u.). Tails of the Bragg peaks are picked up 
at these positions. In addition, diffuse intensities extend- 
ing out from the {100} Bragg peaks along the [110] and 
[110] directions are clearly observed. Around the four 
{110} peaks ((110), (110), (110), and (110)), we found 
that the diffuse scattering intensity is extending more in 
the transverse direction. For example, at the (110) peak, 
the diffuse scattering is strong along the [110] direction. 
These (110) type diffuse intensity distribution is in good 
agreement with previous x-rayi^i^ and neutror&ifi dif- 
fuse measurements. 

When moving to larger Q positions, the Ewald sphere 
starts to move further away from the (HKO) plane in 
the reciprocal space, and the results reveal more inter- 
esting behavior. Around the (200) position, SL ~ 0.09 
r.l.u., two spots are observed, split in the longitudinal 
direction along (100). Taking into account the curvature 
of the Ewald sphere, the positions of the two spots can 
be calculated to be (2.09,0,-0.09) and (1.91,0,-0.09). One 
can see that (2.09,0,-0.09) = 0.09x (1,0,-1) + (2,0,0), and 
(1.91,0,-0.09) = -0.09x (1,0,1) + (2,0,0). It is possible 
that those could be diffuse intensities extending out from 
the (200) peak in the [101] and [101] directions. In addi- 
tion, the [110] and [110] type diffuse intensities can still 
be vaguely seen around the (200) peak. Although these 
diffuse intensities are much weaker than the Bragg peak 
intensity, they are also much broader in reciprocal space. 
This is why when the Ewald sphere is further away from 
the Bragg peak at (200) that that at (100), the sharp 
Bragg peak itself is not observed, but small traces of the 
[110] and [110] type diffuse intensities still remain. 

Around the (210) position, in addition to the two spots 
split along (100) direction, we can see vaguely two weak 
spots split along the (010) direction too. This can be bet- 
ter seen at the (220) position, where the Ewald sphere 
is further away from the (HKO) plane, with SL ~ 0.18 
r.l.u. Four spots are observed, positioned at around 
(2.18,2,-0.18), (1.82,2,-0.18), (2,2.18,-0.18), and (2,1.82,- 
0.18). Similar to the two spots around (200) position, 
these four spots are possibly intensity rods extending out 
from the (220) peak, along the [101], [101], [011], and 
[011] directions, and cutting through the Ewald sphere 
with a none-zero L value. 

In order to verify this, more measurements were per- 
formed by tilting the sample, so that diffuse scattering in- 
tensities at different SL values away from the Bragg peaks 
can be studied. In Fig.0J images taken at T = 15 K and 
500 K, as well as different sample tilt angles arc shown. 
At T = 15 K, sample tilt 9 = 0°, the result (Fig. H 
(a)) looks very similar to that at T = 200 K, shown in 
Fig. When the sample was tilted by 2.6° (Fig. H (b)), 
the (200) Bragg peak was actually right on the Ewald 
sphere, and the CCD was saturated at that position. The 
vertical lines in the image are artifacts of the CCD sat- 
uration. The (200) peak, however, was moved further 
away from the Ewald sphere. The sphere is therefore 
cutting at a larger SL ~ 0.18 around the (200) posi- 
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FIG. 3: CCD image taken from PZN-8PT, at T=200 K. The incident x-ray beam is along the [001] direction. The measurements 
are diffractions coming from very close to the reciprocal lattice (HKO) plane. 



tion. The two spots around (200) on the image clearly 
split further apart than that measured without sample 
tilt (Fig. E|(a)). They are roughly positioned at (2.18,0,- 
0.18) and (1.82,0,-0.18). This confirms that these two 
spots are indeed results from Ewald sphere cutting the 
intensity "rods" coming out from the (200) Bragg peak 
along the [101] and [101] directions. The sample tilt also 
moved the (100), (110) and (110) peaks (on the top part 
of the image) further away from the Ewald sphere. The 



fine features around these positions are also consistent 
with those around (200) and (220). Two longitudinally 
split spots at (100), and four spots at (110) and (110). In 
Fig- E] (a), these Bragg points are simply too close to the 
Ewald sphere for any fine features to be observed. The 
(220) and (220) peaks (on the bottom part of the im- 
age) are, however, moved closer to the Ewald sphere in 
Fig- El (b) then in (a). Here we see the four spots moving 
closer in, also in good agreement with what one would 



■5 



0=0 8=2.6 




PZN-8PT, (HKO) zone 

FIG. 4: CCD image taken from PZN-8PT, at T=15 K and 500 K. The sample tilts are 9 = 0° and 2.6°, respectively. The 
incident x-ray beam is along the [001] direction. The measurements are diffractions coming from very close to the reciprocal 
lattice (HKO) plane. 



expect from the Ewald sphere cutting intensity "rods" 
along [101], [101], [011], and [011] directions. 

One important fact to note is that around the (200) 
position, no intensity rods along the [011] and [011] was 
observed. Similarly, around the (020) position, no inten- 
sity rods along the [101] and [101] directions are present. 
In x-ray and neutron scattering measurements, the dif- 
fuse scattering intensity resulting from correlated polar- 



izations (atomic shifts) is proportional to |Q • e\ , where 
e is the polarization vector. It is therefore evidential that 
the polarizations contributing to the [011] and [011] type 
diffuse scatterings are perpendicular to Q =(1,0,0). 



e?=o" 



0=3,6" 






(100) 

— Jwi 

PZN-8PT, T=15K 
(HKK) zone 

FIG. 5: CCD image taken from PZN-8PT, at T=15 K. The sample tilts are 6 = 0°, 2.6°, and 5.2°, respectively. The incident 
x-ray beam is along the [Oil] direction. The measurements are diffractions coming from very close to the reciprocal lattice 
(HKK) plane. 



B. (HKK) zone 

Diffuse scattering measurements were also performed 
in a different zone, the (HKK) zone. In Fig. results 
are shown at T = 15 K, with the sample rotated by 45° 
around the vertical axis ([100] direction). Now the x-ray 
beam is incident in the [Oil] direction, and the Ewald 
sphere is almost parallel to the reciprocal (HKK) plane, 



i.e., the plane defined by the [100] and [011] vectors. 

With sample tilt 9 = 0°, we see four spots around the 
(200) position as expected, where the Ewald sphere is 
about 5q =(0,-0.064,0.064) away from the (HKK) plane. 
These four spots are roughly positioned at (2.13,0,0.13), 
(2.13,-0.13,0), (1.87,0,0.13), and (1.87,-0.13,0). It is en- 
tirely consistent with the results from the (HK0) zone 
measurements, i.e., intensity rods extending out in the 
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[101], [101], [110], and [110] directions from the (200) 
peak, being cut by the Ewald sphere. With tilt 9 = 2.6° 
(see Fig. 03(b)), these four spots move further out, as the 
Ewald sphere is moving further away from the (HKK) 
plane around (200). No central peak is seen at the cen- 
ter of those four spots (Q —(2. ,-0.64, 0.64)), confirming 
the absence of diffuse intensity rods in the [Oil] direc- 
tion around the (200) position. In addition, one starts 
to see similar features around the (100) position as well. 
With further tilting, 6 = 5.2° (see Fig. EI (c)), the intensi- 
ties around (200) has become very weak, and the further 
parted four spots can hardly be observed. On the bot- 
tom of the image, four spots also appear around the (200) 
position. There the Ewald sphere actually went passed 
the (HKK) plane, with Sq =(0,0.064,-0.064), so similar 
patterns are seen there compared to Fig. El (a). 

Diffuse patterns around other Bragg positions are all 
consistent with this (110) type diffuse rods. For exam- 
ple, if we look at around the (220) position, a strong peak 
appears there. The Bragg structure factor of PZN (220) 
and (200) are the same order of magnitude, yet the Ewald 
sphere is further away from the Bragg peak at (220) than 
(200) (because of the larger Q). Since the (200) Bragg 
is not observed, neither should the (220) Bragg. In fact, 
that "peak" is not the (220) Bragg peak, but instead the 
diffuse scattering intensity coming out from the (220) 
peak along the [011] direction (perpendicular to the im- 
age plane). 



C. 3-D diffuse intensity distribution 

Our results can establish unambiguously the geometry 
of the diffuse scattering intensities distributed in the re- 
ciprocal space. Based on all the information obtained, 
(i) diffuse intensities are (110) type rods around Bragg 
peaks; (ii) the intensities of different (110) diffuse rods 
can vary between different Bragg reflections; (iii) at some 
Bragg reflections, certain (110) intensity rods are strong, 
and certain (110) rods may even be absent, the 3-D distri- 
bution of the diffuse scattering intensity can be sketched 
in Fig. El The 3-D sketch is shown in the way that the 
(001) plane ((HK0) plane) is parallel to the paper. Here 
we plot the equi-intensity surface of the diffuse scattering 
around various Bragg peaks. The contrast and sharpness 
of the diffuse intensity rods have been exaggerated in or- 
der to illustrate the geometry in a simple and effective 
way. 

In a quantitative way, it describes our results very 
well. For example, around the (100) peak, the [110], 
[110], [101], and [101] diffuse rods are of equal inten- 
sities, but the [011] and [011] rods are absent; around 
the (110) peak, the [110] diffuse rod is much more in- 
tense than the [101], [101], [011], and [011] rods, but 
the [110] rod is absent. A cross section of the 3-D dif- 
fuse scattering distribution in the (HK0) plane gives the 
"butter-fly" shaped pattern around the (100) peak, and 
the transverse [110] type diffuse around the (110) peak, in 




FIG. 6: Sketch of the diffuse scattering distribution in the 
3-D reciprocal space around (100), (110), (111), (010), and 
(011) reciprocal lattice points. 

good agreements with previous neutron scattering mea- 
surements&iS. More detailed model calculations will be 
presented in the next section. 

IV. DISCUSSIONS 
A. Static or dynamic? 

One disadvantage of x-ray diffraction measurements is 
the coarse energy resolution, so that many of the low 
energy phonon modes are also included together with 
the elastic component. Some examples of phonon con- 
tributions dominating the x-ray diffuse intensity distri- 
bution can be found in similar ferroelectric perovskite 
compounds BaTiOj 22 and KTaOs 23 . A very important 
question to resolve is whether these intensity rods ob- 
served in PZN-xPT diffuse are coming from phonon con- 
tributions (thermal diffuse) or static distortions (static 
diffuse) . 

Recently, neutron diffuse scattering measure- 
meru - a 8 i io i i3 4 24 nave Deen performed on the relaxor 
systems with good energy resolutions so that phonon 
intensities can be easily separated out. These results 
have directly confirmed that the (110) type diffuse 
intensities in the (HK0) plane are elastic. 

The temperature dependence of the diffuse scattering 
intensity is also consistent with the static nature. In 
Fig. the temperature dependence of the diffuse scat- 
tering intensity at (2.09,0,-0.09) (the top one of the two 
spots observed around the (200) position) are plotted for 
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PZN-xPT Q=(2.09,0,-0.09) 



to 
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FIG. 7: Temperature dependence of the diffuse scattering 
measured at (2.09,0,-0.09) for PZN, PZN-4.5PT and PZN- 
8PT. 



PZN, 4.5PT, and 8PT. At T below T c , the intensity de- 
creases slowly and monotonically with increasing T for 
all three compounds. In addition, the diffuse intensity 
does not reach zero at the ferroelectric phase transition, 
but rather at some temperatures above Tc- All these are 
in agreement with the elastic diffuse scatterings expected 
from PNR in relaxor ferroelectric systems. 

Phonon scattering intensities, on the other hand, are 
proportional to ksT / '(fnu) 2 , when taking into account of 
the Bose factor ksT/hu and the 1/ui factor. In PZN 9 
and PMN 25 , the zone center TO phonon energy actually 
increases with cooling for T < Tc- Measurements on TA 
phonons at small q in the same systems 9,2 — also confirmed 
that the TA phonon modes behave normally below Tc, 
i.e., no significant broadening or softening. Therefore, 
both TO and TA phonon contribution to the x-ray diffuse 
scattering should increase with increasing T, at least for 
the low temperature range T < Tc- This is in contrast to 
the data shown in Fig. [7| As a comparison, x-ray diffuse 
scattering measurements with similar set-ups on single 
crystals of Si have been performed by us, and by Holt et 
alm^. There the phonon contribution dominates, and the 
thermal diffuse intensity clearly increases with heating. 

Another clue is provided by comparing the structure 
factors calculated based on previous results. Previous 
measurements by You et al. obtained similar results near 
the (300) Bragg peak. It was attributed to contributions 
from ferroelectric soft transverse optic phonons propa- 
gating along the (110) directions, and thus concluded to 
be only strong around the odd integer Bragg reflections 
(H + K + L = odd). In our measurements, however, 



the diffuse scattering intensities around the (200), (110), 
and (220) Bragg peaks are reasonably strong. The static 
diffuse scattering, acoustic and and soft optic phonon 
structure factors at different Bragg reflections are shown 
in Table. [I] Apparently the diffuse scattering intensities 
observed in our measurements are more consistent with 
the static diffuse structure factors, which do not vary 
much between different Bragg peaks in our measurement 
range. It is also qualitatively in agreement with measure- 
ments by Takesue et alM on PMN in the (HK0) plane. 
The soft phonon structure factors, however, vary a lot 
and are very strong at (210) and (300). Measurements 
around those peaks did not show any particularly strong 
diffuse intensities (see Figs. |3]and QJ. 

These facts strongly suggest that for temperatures T < 
Tc, the main contribution to our x-ray diffuse scattering 
is coming from static lattice distortions, namely, those 
from the polarized PNR, instead of thermally activated 
inelastic phonon scatterings. We can therefore study the 
low temperature diffuse scattering distributions to bet- 
ter understand the static polarizations (atomic shifts) 
and correlations in the PNR. With increasing temper- 
ature, the phonon contributions will eventually increase 
and dominate. In fact, at 700 K, with the disappearance 
of the two spots around (200) position, we started to pick 
up intensities around (2,0,-0.09), which is very likely com- 
ing from thermally activated phonons propagating in the 
[001] direction. 

In the next two subsections, we will start at provid- 
ing a simple phenomenological model describing the real 
space correlation of polarizations, that gives rise to the 
(110) rod-type diffuse intensities. Based on this model, 
the shape of diffuse scattering in these compounds can be 
well calculated within good agreement to the experimen- 
tal results. We will then suggest one possible origin of 
these elastic diffuse scatterings from the aspect of lattice 
dynamics - condensation of soft optic phonons into static 
distortions. 



B. Planar correlations - "pancake" model 

Having determined the static nature of these observed 
(110) type diffuse intensity rods, we now consider a 
simple model that can reproduce such diffuse intensi- 
ties. Generally, a rod type structure in the reciprocal 
space corresponds to a planar structure in the real space. 
Therefore, the (110) type diffuse scattering intensity can 
be a result of polarizations correlated in the {110} planes 
in the real space. There are totally six {110} planes, cor- 
responds to the six (110) diffuse intensity rods. Further- 
more, the polarizations should have the same periodicity 
as the lattice itself, because all these diffuse scattering in- 
tensities peak at the Bragg positions, instead of forming 
some super-lattice peaks. 

The problem remaining is to determine the directions 
of these polarizations with planar type correlations. The 
polarizations can be either in-plane, or out-of-plane, as 
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TABLE I: X-ray scattering structure factors for PZN. The diffuse scattering structure factors are based on the atomic shift 
values derived from measurements by Vakhrushev et alw— on PMN. The soft optic phonon structure factors are calculated 
based on a ratio of S = 1.5 between the Last Mode and Slater mode contributions, as determined by Hirota et alt~ 
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FIG. 8: Polarizations confined in the "pancake" shaped 
spaces. The correlation in the plane (diameter) is much larger 
than the correlation perpendicular to the plane (thickness), 
(a) and (b) are two different models where the polarizations 
are in the plane, or out of the plane, respectively. 



shown in Fig. [S] Here the polarizations are confined in 
the "pancake" shaped spaces, indicating the planar type 
correlations. One can easily find that out-of-plane type 
polarizations as shown in Fig. [81(b) do not yield the type 
of diffuse intensity pattern around the (110) peak. There 
when measured in the (HK0) plane, the diffuse intensity 
is mostly transverse along the [110] direction, and the 
[110] type intensity rod is absent. Since Idiff oc |Q • e| 2 , 
polarizations e correlated in the (110) plane, that give 
rise to the [110] type diffuse, must be perpendicular to 
Q =(1,1,0) for this intensity to be zero. Only in-plane 
type polarizations can satisfy this condition. 

Nevertheless, the polarizations still can not rotate 



O 



FIG. 9: Real space model indicating the correlation of the 
polarizations (polarizations) in the system. Here the [110] 
polarizations are correlated in the (110) plane. 



freely in the {110} planes. The exact direction of these 
polarizations can be determined by studying the diffuse 
intensity around the (200) peak. Here the diffuse scatter- 
ings along [011] and [Oil] directions are absent. So, for 
the [011] type diffuse, e must be perpendicular to both 
Q =(2,0,0) (so that Q • e = 0) and the out-of-plane vec- 
tor [011] (so that e is in-plane). This immediately leads 
to e =[011] . Similar analysis can be used to derive the po- 
larizations directions associated to all the six (110) type 
diffuse scatterings. 

In summary, we find that the real space structure cor- 
responding to the (110) rod- type diffuse intensity is {110} 
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planar correlations of in-plane (110) polarizations. More 
specifically, there are polarizations along the [110] direc- 
tion, correlated in the (110) plane; polarizations along 
the [110] direction, correlated in the (110) plane; polar- 
izations along the [101] direction, correlated in the (101) 
plane; polarizations along the [101] direction, correlated 
in the (101) plane; polarizations along the [011] direction, 
correlated in the (011) plane; and polarizations along the 
[011] direction, correlated in the (011) plane. Fig. [§|pro- 
vides a schematic illustration of one of these six planar 
correlations in real space. The x-ray diffraction measure- 
ments are less sensitive to the lighter atoms such as O, 
so we only illustrated possible motions of the heavier A 
(Pb 2+ ) and B (Zn 2+ /Nb 5+ /Ti 4+ ) site atoms. The mo- 
tions of the A site and B site can be parallel (acoustic) 
or antiparallcl (optic) , or a combination of the two. 




FIG. 10: Model calculation of the diffuse scattering distribu- 
tion in the 3-D reciprocal space around (100), (110), (111), 
(010), and (011) reciprocal lattice points. 

The diffuse scattering intensity can be therefore cal- 
culated based on this simple model. At Q = G + q, 

7 rfi// (Q) = A|F dl// (G)| 2 £ |Q • £i | 2 ^X_. 

(1) 

IdiffiQ) is the sum of contributions from the six {110} 
planar type correlations, where |-Fdi//(G)| 2 is the diffuse 
scattering structure factor at the Bragg reflection G; £j 
is the polarization vector correlated in the plane (one of 
the six (110)); T\\ = and Tj_ = l/£±, are the in- 

verse of the correlation length in- and out-of-plane; 
and q±i are the in- and out-of-plane components of q. 
Here we used a product of Lorentzian functions because 



an exponentially decaying correlation in real space, both 
in- and out-of-plane, has been assumed. For a square- 
function cut-off type correlation caused by a finite "pan- 
cake" shape, with thickness described by £j_ and diameter 
described by the reciprocal space intensity distribu- 
tion will take the functional forms of Gaussian functions 
instead. However, the shape of the diffuse scattering in- 
tensity will not change qualitatively. 

In the simulation, the in-plane correlation length was 
assumed as 20 lattices, 4 times larger than the out-of- 
plane correlation length (5 lattices - estimated from re- 
sults of neutron diffuse measurements on PMN by Xu et 
al&). The contributions from the six {110} planar corre- 
lations are weighed by the |Q • e| 2 factor and summed up. 
The result is shown in Fig. 1101 Here the equi-intensity 
surface of the diffuse scattering intensity is plotted. Note 
that for best visual effects, the static diffuse structure 
factors in Tableland the |Q| 2 factor were taken out in 
constructing Fig. ^3 m order to illustrate the diffuse in- 
tensity distributions around different Bragg peaks with 
the same quality in the same figure. In the exaggerated 
version of this model, as shown in Fig.|S| the in-plane cor- 
relation length was chosen to be 20 times larger than the 
out-of-plane one to produce those "sharp" diffuse rods. 

In Fig. the simulated diffuse scattering intensities 
around the (100), (010), (110), (111), and (011) Bragg 
peaks are plotted. When measured in the (HK0) plane, 
the cross section of the 3-D diffuse scattering intensi- 
ties gives the "butter-fly" diffuse pattern around (100) 
and (010) peaks, as well as the transverse [110] type dif- 
fuse around the (110) peak. By moving away from the 
(HK0) plane, with a small 8L, one starts to see the two 
spots split in the longitudinal direction around the (100) 
Bragg peak due to the [101] and [101] type diffuse inten- 
sities; and four spots around the (110) Bragg peak due 
to the [101], [101], [011], and [011] type diffuse intensi- 
ties. A simulation of diffuse intensities in the (HKSL) 
plane, at different SL values is shown in Fig. ^] It is 
in good agreement with our measurements in the (HK0) 
zone. Similarly, one can also compare the simulation in 
the (HKK) zone to our measurements. We found that 
calculations based on this model can explain all of our 
results, as well as almost all of the previous results from 
x-rajii^i^ and neutron&iflii^ diffuse scattering measure- 
ments on PMN and PZN-xPT. Although it does not yet 
provide a clear physical picture on the origin of these cor- 
related polarizations, it is very useful in calculating and 
predicting diffuse scattering distributions in these relaxor 
compounds. 

One interesting problem is that all of the three com- 
pounds in our measurements have rhombohedral type po- 
larizations below Tc- Recent NMR measurements^ on 
PMN also confirms the (111) type polarization in the sys- 
tem. How do the (111) type polarization and the corre- 
lated (110) type polarization coexist in the same system? 
In fact, our model is not in contradiction to the average 
(111) type polarization. Alternatively, the (111) type 
polarization can be decomposed into three (110) com- 
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FIG. 11: Model calculation of the diffuse scattering distri- 
bution in the (H,K,0), (H,K,0.05) and (H,K,0.1) reciprocal 
planes (i.e. (001) planes, at different small L values), around 
the (100) and (110) reciprocal lattice points. 

ponents. The average (111) type polarization can be a 
result of global averaging of (110) type polarized "pan- 
cake" entities. Another possible scenario could be the 
following: when three different {110} type "pancakes" 
(e.g., (110), (011), and (101)) cross each other, the region 
shared by those three "pancakes" would have a combi- 
nation of three (110) polarizations - resulting in a (111) 
type polarization. This shared region could be a (111) 
polarized PNR, and the size of this region is then defined 
by the thickness of the three "pancakes", i.e., the out-of- 
plane correlation length, which is in the order of one to 
few nanometers (a few lattice units). 

C. Condensation of (110) polarized soft optic 
phonon modes 

Ferroelectric polarizations are optic type distortions, 
associated with different atoms moving opposite to each 
other in the unit cell. It is well accepted that the PNR 
are formed from the condensation of the ferroelectric soft 
optic phonon mode. Previous measurements show that 
the transverse optic (TO) phonons become overdamped 
for Tc < T < Td for q smaller than a certain value, 
called the "waterfall" wave-vector q w j22£!^. The "water- 
fall" effect was observed for q along both [001] and [011] 



directions^, and is believed to be intimately related to 
some length scales defined by the PNR. We can therefore 
suggest one possible origin of the rod- type diffuse scatter- 
ing intensities. If the (110) polarized optic phonon modes 
are particularly soft when propagating along the perpen- 
dicular (110) directions, they can then condense into the 
(110) type elastic diffuse intensities upon cooling. This 
occurs around T = Td, when the soft optic mode becomes 
overdamped and start to condense into static distortions. 
The shape of these static entities are then determined by 
the energy dispersion, i.e., phonon velocity, of these (110) 
polarized modes. Take the [110] polarized mode for ex- 
ample, it is likely to be softer along the [110] direction, 
i.e., the optic phonon propagates with a smaller velocity 
along the [110] direction; and harder in the (110) plane. 
When this mode becomes overdamped at Td, it then con- 
denses accordingly, with a shorter length scale along the 
[110] direction than in the (110) plane, resulting in the 
(110) planar correlation of the [110] polarizations, and 
the "pancake" shaped entity. 



Since PNR are results of soft TO phonon condensa- 
tion, why do the static diffuse scattering and soft optic 
phonons have different structure factors? This discrep- 
ancy has been very well explained by Hirota et al. with 
the phase-shifted PNR— model. They re-examined the 
atomic shift values derived from the neutron diffuse scat- 
tering measurements by Vakhrushev et al. 2 " 7 on PMN, 
and found out that those values do not satisfy the center 
of mass condition. The mass center of the atoms in the 
PNR are displaced from the surrounding lattice along 
the polarization direction. This displacement is called 
the "uniform-phase-shift" , which is an acoustic type of 
atomic motion. The ferroelectric soft phonon mode only 
condense into the center-of-mass portion of the atomic 
shifts in the PNR, and therefore have different structure 
factors than the overall static atomic shifts in the PNR. 



In our measurements, we were able to probe dif- 
fuse scattering intensities distributed in three dimensions 
around many Bragg peaks simultaneously from PZN- 
xPT single crystals. The shape of the diffuse scattering 
in the reciprocal space can be described as (110) "rod" 
type intensities, resulting from {110} planar correlations 
of (110) polarizations in the real space, with the in-plane 
correlation length about a few to ten times larger than 
the out-of-plane correlation length. On the other hand, 
limited by the measurement resolution, we can not accu- 
rately determine the values of correlation length in and 
out of the correlation planes. In addition, we also do 
not have enough information to demonstrate whether the 
in-plane correlations are isotropic, or have certain pref- 
erences. Further studies on these relaxor systems are 
required to complete these details and to provide a clear 
physical picture to explain the origin of the planar cor- 
relations. 
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